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We present the final results of a search for lightly ionizing particles using the entire cosmic ray data set 
the MACRO detector collected during its 1995-2000 run. Like the original search performed with the data of 
1995-96, this search was sensitive to fractionally charged particles with an electric charge q as low as | and 
with velocities between approximately 0.25c and c. The efficiency of this search was ~ 70 % for q = | and 
increased rapidly to 100 % for higher charges. No candidate events were observed. This corresponds to a 90 % 
C. L. upper limit on their isotropic flux of 6.1 x 10~ 16 cm' 2 sec -1 sr _1 which represents the most stringent 
experimental limit ever obtained. 
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PACS numbers: 14.80.-j, 96.40.-z 



I. INTRODUCTION 

Within the Standard Model of electroweak and strong in- 
teractions electric charge quantization remains unexplained. 
Since Millikan's oil drop experiment demonstrated electric 
charges in Nature come in discrete units, numerous experi- 
menters have refined the original determination of the unit 
charge (e) and searched for free particles with fractional 
charge. This interest was intensified with the proposal of the 
quark model by Gell-Mann Ql and Zweig |y] in 1964 and 
its undisputed success in explainingthe deep inelastic scatter- 
ing experiments in the late 1960's y|]. This model -now part 
of the Standard Model- does require that single quarks come 
with a fractional charge, however, it does not allow them to 
exist as free particles. Quarks are confined in color-neutral 
baryons and mesons that carry integer electric charges. 

Searches for fractionally charged particles have thus been 
carried out over the last decades in bulk matter, at acceler- 
ators and in cosmic rays IH 0. S 01 without any ev- 
idence for their existence. The observation of such parti- 
cles would be direct evidence of physics beyond the Stan- 
dard Model. Fractionally charged particles are easily acco- 
modated in Grand Unified Theories (GUT) of the electroweak 
and strong interactions. Within this framework the quantiza- 
tion of the electric charge is also explained as a consequence 
of the non trivial commutation relations between the opera- 
tors in the theory. Simple extensions of the SU(5) unifica- 
tion group allow for | and ^ char ges fiol [Till while larger 
groups such as SU(8) d, SO(14) (H and SO(18)Jp al- 
low for | ones. Many popular superstring models 11511 yield 
stable particles with charges down to | or even smaller. Fi- 
nally, some spontaneously broken QCD theories predict the 
possibility of primordial unconfined quarks and gluons fl^l 
contained in superheavy quark-nucleon complexes with large 
noninteger charge. 

Previous searches for fractionally charged particles in the 
penetrating cosmic radiation include the search using the first 
year of MACRO'S running fl7ll as well as searches with 
the Kamiokande-II fl8il and LSD fl9il experiments. Un- 
til this search the best 90 % C. L. upper flux limits were 
obtained by the Kamiokande experiment and are 2.1 and 
2.3 x 10~ 15 cm -2 sec -1 sr _1 for charges | and respec- 
tively. Here we present the final results of the search for par- 
ticles with fractional charge from | to y in the penetrating 
cosmic radiation using the entire dataset of the MACRO ex- 
periment. The search method is based on the assumption of 
a reduced ionization and atomic excitation energy loss rate of 
a fractionally charged particle with respect to that of a unit 
charge. Since the energy loss rate of a particle of charge Q 
is proportional to Q 2 , a fractionally charged particle is ex- 
pected to lose energy, by excitation and ionization, at a much 
lower rate than a minimum ionizing particle of the same ve- 
locity. The fractionally charged particles are therefore also 
called Lightly Ionizing Particles (LIPs). This search would 
not have been sensitive to fractionally charged particles that 



may interact also strongly (e.g., standard model quarks) as 
these particles would not be able to penetrate large amounts 
of material. 



II. EXPERIMENTAL SETUP 

The MACRO experiment, described in detail in I20I ElTl . 
was a multipurpose underground detector located in Hall B 
of the Gran Sasso National Laboratory of the Italian Insti- 
tute of Nuclear Physics (INFN). It was optimized for the 
search for magnetic monopoles 1221 l23ll but it also presented 
excellent capabilities for studying atmospheric neutrino os- 
cillations 1241 l25l | 2 ja|. cosmic rays 1271 12811 . astrophysi- 
cal point sources neutrinos from gravitational col- 

lapses ll3ll l32ll . The detector had a modular structure; it 
was composed of six supermodules each of which was di- 
vided into a lower and an upper part ("attico"). All super- 
modules were equipped with three detector sub-systems: lim- 
ited streamer tubes, liquid scintillator counters and nuclear 
track detectors. The overall dimensions of the apparatus were 
76. 6x12x9. 3 m 3 and the acceptance for an isotropic flux of 
particles was fs 10 4 m 2 sr. The detector was active in various 
different configurations from the fall of 1989 until December 
2000. 

Lightly ionizing tracks in MACRO would be signatures 
for fractionally charged particles. They could be identified 
both in the streamer tube and the scintillator system. The 
streamer tube system provided a three-dimensional recon- 
struction of the particle's trajectory. This system was more 
than 99% efficient in forming a track upon the passage of a 
fractional charge |e or greater The scintillator sys- 

tem provided redundant measurements of the particle's energy 
loss. Time-of-flight information with sub-nanosecond accu- 
racy was also recorded. A Lightly Ionizing Particle (LIP) 
trigger was formed in hardware by combining MACRO'S pri- 
mary streamer tube trigger I17ll2lll34ll with a coincidence of 
MACRO'S lowest e nerg y threshold scintillator trigger, called 
PHRASE 1201 12 ll Bill , coming from at least three detector 
faces. The coincidence window of 400 ns among the scintil- 
lator face triggers was setting the minimum particle velocity 
the trigger hardware was sensitive to. This was at about 0.25c 
for most of the detector's acceptance. PHRASE provided a 
trigger based on a single scintillator counter with a threshold 
at approximately 1.2 MeV. Given that the Landau peak for 
cosmic ray muons was approximately 38 McV for a typical 
20 cm pathlength within a scintillator counter (see fig. ^ the 
1.2 McV energy threshold set a minimum fractional charge 
of \e that MACRO could detect. A typical plot of the mea- 
sured trigger efficiency of the low-energy PHRASE trigger as 
a function of energy is shown in fig. |2] as obtained with the 
1995-96 data 0. The low energy threshold of the PHRASE 
was periodically checked throughout the five-year run using 
natural radioactivity decay lines and it was maintained con- 
sistently at the 1 — 1.2 McV level during the entire run. 
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FIG. 1: Typical energy loss rate in a MACRO liquid scintillator 
counter. The superimposed fit is a Landau distribution folded with 
the photoelectron fluctuations. One may notice the low dE/dx tail of 
the energy loss distribution that extends within the LIP signal region. 



The LIP circuitry provided the hit pattern information of 
the scintillator counters that exceeded the 1.2 MeV threshold. 
For all counters involved in a LIP trigger their photomultiplier 
(PMT) wave forms recorded on MACRO'S 200 MHz wave 
form digitizer (WFD) system were read. An analysis of these 
wave forms could then provide the localization of the event 
along the scintillator counter as well as the reconstruction of 
the energy released. Both of these when combined with the 
three-dimensional track information provided by the streamer 
tubes could effectively reject any radioactivity or cosmic ray 
muon signal at a level of better than 1/10 6 leaving a very low 
background search for LIPs. This was the approach in our first 
search for LIPS using MACRO'S 1995-96 data O. 



m. DATA ANALYSIS 

The MACRO scintillator subdetector was equipped with 
three independent systems able to record the time of flight, 
position and energy loss of a lightly ionizing particle as it 
crossed the apparatus: the primary muon trigger (called ERP), 
the primary gravitational collapse trigger (called PHRASE) 
and the Wave-Form-Digitizer system (WFD). While all three 
had comparable resolutions 1171 l20ll2lt I2lll25l 1321 for mini- 
mum ionizing particles their different thresholds implied sen- 
sitivity to as low as ~ 4p for the ERP (given its ~ 15 MeV 
threshold), ~ ^ for the PHRASE (given its ~ 7 MeV thresh- 
old), and ~ | for the LIP-triggered WFD system. The LIP 
trigger was designed to trigger on cosmic ray muons in addi- 
tion to lightly ionizing particles. We have thus used the ERP 
and PHRASE systems in order to identify and reject the cos- 
mic ray muon component that was present in our LIP trig- 
ger data set. This allowed the wave form digitizer system to 
be invoked only when the ERP and PHRASE had failed to 
record a candidate event presumably because of their energy 
thresholds and triggering inefficiencies relative to LIP/WFD 
system. The energy thresholds of the PHRASE and ERP sys- 
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FIG. 2: Measured trigger efficiency of the low-energy PHRASE 
trigger (i.e., LIP trigger) as a function of the energy released in the 
liquid scintillator counters. Some measured efficiencies were larger 
than 100 % because the normalization factor used was an estimate of 
the true normalization as a function of energy. 

tems were continuously monitored using natural radioactivity 
and cosmic ray muons. They were periodically adjusted in 
order to account for any electronic drifts or photomultiplier 
tube gains. This guaranteed their consistency and maintained 
our ability to trigger and reconstruct tracks down to | charges 
throughtout the entire data taking period. 

The data set for this search comes from the 5-year run of 
the MACRO detector in its final configuration from July 1995 
to December 2000. Run quality criteria were first applied; 
they required the MACRO detector to be running in its full 
configuration (i.e., all six supermodules, excluding counters 
with calibration errors) and without any serious acquisition 
problems (e.g., high dead time, repeated hardware errors, high 
voltage power supply glitches or high voltage errors what-so- 
ever). The integrated live-time was 1320 days and during that 
period 18.3 x 10 6 LIP triggers were collected. The analysis 
of these triggers proceeded as follows: 

1) Using the streamer tube hit information, we first re- 
quired a single streamer tube track reconstructed. In 
addition, using the LIP scintillator counter hit informa- 
tion, we required hits to be present in no more than four 
scintillator faces and within six scintillator counters in 
the same face. Given the geometry of the detector, no 
single candidate LIP track could result in such a scintil- 
lator counter hit pattern. This cut primarily rejected cos- 
mic ray muons accompanied by electromagnetic show- 
ers. 

2) Using the streamer tube track parameters we recon- 
structed the event longitudinal position along the scin- 
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tillator counters as well as the path length of the cross- 
ing particles in them. We required the path length 
to be between 13 cm and 70 cm and the hit position 
along the counter within the central 10.8 m part of it 
(i.e., rejecting tracks within the final 10 cm of a scin- 
tillator counter). Both these cuts eliminated track ge- 
ometries that were more susceptible to energy recon- 
struction errors either due to low photoelectron statis- 
tics (short pathlengths) or to poorly calibrated PMT re- 
sponses (near-PMT geometric response). When avail- 
able, we required that the position along the counter 
provided by the scintillator timing was in agreement 
within 80 cm (6 a ) with the position reconstructed by 
the streamer tube system. Following these geometrical 
cuts, the detector acceptance was about 3300 m 2 sr for 
an isotropic flux of particles. 

3) For every LIP counter that was intercepted by a track 
that fulfilled all cuts so far, we assigned the energy loss 
as calculated for the same scintillator counter(s) by the 
ERP muon system in the absence of which this informa- 
tion was sought in the PHRASE gravitational collapse 
system. In lack of any energy loss information from 
these two systems a zero energy loss was assigned to 
that hit; this makes a priori the best LIP candidate event 
in MACRO. 

This last stage of the LIP analysis enabled us to perform 
a measurement of the particle energy loss rate dE/dx for all 
cases that the LIP event was accompanied by a muon and/or 
the gravitational collapse trigger. The maximum dE/dx 
among the counters involved is plotted in figure [3] All tracks 
in this analysis were accompanied by a muon and/or gravita- 
tional collapse hit. This allowed us to establish the particle's 
dE/dx without making use of the wave form digitizer sys- 
tem. Since energy information from both the PHRASE and 
ERP systems was used in this search, we verified that in the 
relevant energy range (104- 100 McV) the two systems agreed 
to within 20 % or better in more than 95 % of the liquid scin- 
tillator counters (see fig.EJ. 

IV. RESULTS AND DISCUSSION 

The expected signal region for this LIP search was set be- 
low 1.1 MeV/cm for the maximum energy loss rate -when 
measured by the ERP- in any of the scintillator counters in- 
tercepted by a streamer tube track. Since the energy loss 
was measured by PHRASE whenever ERP information was 
not available and in order to account for the potential 20% 
mismatch in their energy estimates, the signal region was ex- 
tended up to 1.35 MeV/cm whenever a measurement was per- 
formed by PHRASE. Above this level the cosmic ray muon 
energy loss spectrum does not permit any identification of 
fractionally charged particles thus setting an upper limit of 
sensitivity to approximately ^ . 

As one can see in fig. [3] there is one event (run 15871, 
event 5649) that appears in the signal region. It corresponds 
to a maximum energy loss of 0.66 McV/cm, i.e., about 20 % 
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FIG. 3: Energy loss as measured by PHRASE for the LIP events 
that passed the track quality and geometry cuts and satisfied the re- 
quirement of a maximum energy loss rate (measured by ERP) less 
than 1.1 MeV/cm. This dataset was comprised primarily of events 
where the ERP system did not trigger at all due to inefficiencies. For 
these events an energy measurement was sought in the PHRASE sys- 
tem. The signal region is in the [0, 1.35] MeV/cm interval. 
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FIG. 4: Comparison between the PHRASE and ERP measured en- 
ergies in the range of 10 -j- 100 MeV. The distribution of the dif- 
ference between the two energies normalized to their average shows 
how well these two agree. The distribution is roughly gaussian with 
a = 0.13. 



lower than what expected for a particle of charge =^ and about 
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a factor of 3 higher than what expected for a particle of charge 
|. Three scintillator counters were involved in this trigger; the 
first in one of the upper vertical layers, the second in the cen- 
tral horizontal layer and the third in the lower horizontal layer. 
There were no ERP triggers involved (suggesting the energy 
released in each of the counters was below 15 MeV) and only 
one PHRASE hit in the vertical layer was present (suggest- 
ing the energy released in the other two counters was below 
7 MeV). The energy loss measured by PHRASE for this hit 
was 13.7 MeV and using the path length in the box provided 
by the tracking we computed a maximum energy loss rate of 
0.66 MeV/cm. The position along the counter for this par- 
ticular box measured by the PHRASE and by the streamer 
tube track geometry were in agreement (within 15 cm). We 
have examined this event by hand relying primarily on the 
wave forms as recorded for all the counters involved in the 
trigger. The apparent amplitude of the recorded wave forms 
was consistent with the energy thresholds for the ERP and the 
PHRASE. Having three scintillator counters involved in the 
trigger we have checked for a consistency in the relative tim- 
ing of them with the crossing of a single particle of constant 
velocity. The relative timing between the counter in the upper 
part of the detector and that in the central part was consistent 
with the passage of a relativistic particle coming from above 
while the relative timing between the box in the lower part of 
the detector and any of the other two hits was consistent with a 
slowly moving upward-going particle. We thus discarded this 
event from the signal region. 

V. CONCLUSIONS 

Over the five years of running of the MACRO detector in 
its final configuration (1995-2000) we observed no candidates 
for LIPs in cosmic ray single tracks. 

Given the detector's acceptance and live time we computed 
an integrated exposure of 3.8 x 10 15 cm 2 ssr which for an 
isotropic flux of particles yields a 90% C.L. upper flux limit of 
6.1 x 10~ 16 cm -2 s _1 sr _1 . The search was fully efficient for 
charges | to ^ . The flux limit rises to higher values for lower 
charges as the detection efficiency drops from 100% to lower 
values. This limit improves the previously published MACRO 
result fl7ll by over an order of magnitude. The previously 
published MACRO search applies to LIPs present in both sin- 
gle and multi-track events while this one (as well as the ones 
from experiments other than MACRO) applies to LIPs present 
in single tracks only. 

This MACRO result can be compared with those obtained 
by the Kamiokande fiHl and the LSD fl^l experiments keep- 
ing in mind that all these are limits on the local flux of LIPs 
at the detector site. Kamiokande and LSD quoted their results 
for two specific values of fractional charge, namely, | and 
while the MACRO result applies to a continuum of charges. It 



should also be noted that all limits were obtained with the as- 
sumption of an isotropic flux of LIPs. This is reasonable if we 
assume that LIPs are produced in the rock around the detector 
(for instance, by cosmic ray muon interactions or some other 
unknown mechanism). However, if LIPs are produced by the 




i i i i i i i i i i i 

6.0 5.0 4.0 3.0 2.0 1.0 

e/q 



FIG. 5: The 90% C.L. MACRO flux upper limits for LIPs 
(solid line) compared with previous limits from MACRO, LSD and 
Kamiokande. 

cosmic ray interactions in the upper atmosphere or impinge 
on the earth from outer space, the assumption of an isotropic 
flux is no longer valid. In this case a detailed physical model 
for the production and propagation of LIPs as well as of the 
response of the detector would be needed in order to derive 
and compare results for the various experiments. In the most 
obvious case, LIPs without enough initial kinetic energy will 
not be able to reach the detector from directions below the 
horizon thus worsening all upper limits by a factor of two. 
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